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ABSTRACT 


This work is focused on the use of biomass waste to feed already existing coal combustion plants as a part of paving 
the way toward the reduction of the environmental impact. The biomass waste supply chain optimization is critical 
to conceive long-term viable projects and deal with the biomass heterogeneous nature and drawbacks to be used 
with coal, i.e. principally high moisture content and low bulk density. This paper studies biomass transportation, 
storage and change of properties (moisture content and hence dry matter, energy density and bulk density) through 
the use of different pre-treatments: (i) torrefaction, (ii) torrefaction combined with pelletization, (iii) pelletization, 
(iv) fast pyrolysis and (v) fast pyrolysis combined with char grinding, which produce a range of very different pre¬ 
treated biomass. The optimization problem is formulated as a mixed integer linear program (MILP) that evaluates 
the net present value and the environmental impact through a life cycle assessment (LCA). The results propose 
location-allocation decision together with the selection/capacity of pre-treatment technologies for each scenario. 
The scenarios contemplate different biomass characteristics, availability and distribution for a supply chain case 
study located in Spain: forest and agricultural woody residues used to replace at least 10% of the total thermal inlet 
power provided by coal in the existing network of thermal plants. 

© 2014 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved. 
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1. Introduction 

The energy sector aims at increasing the use of renewable 
sources and of efficient processes. Countries all over the world 
have become conscious of the biomass potential as a renew¬ 
able source to contribute to the satisfaction of energy needs, 
promote development in rural areas and alleviate climate 
change. Biomass can be of different origins and, among the 
different types, organic wastes are of special concern to avoid 
land use controversy and to propose a safe final disposal 
while generating electricity, fuels and/or heat/cool. Biomass 
can contribute in centralized or large scale, as well as in dis¬ 
tributed or small scale projects, by satisfying different user 
needs. Biomass use in centralized energy systems (CES) needs 


from well-established channels to receive raw material on 
time at the required characteristics, and to distribute the final 
product. This final output can be any bio-product: chemi¬ 
cal compounds, electricity, heat, hydrogen, synthetic natural 
gas and/or liquid fuels (Higman and van der Burgt, 2003). 
As an emerging sector, there is a need to define its basis to 
secure that the biomass processing project as a whole is eco¬ 
nomically, environmentally and socially sustainable with the 
current and future tendencies and policies, from the renew¬ 
able source to the grid (Cherubini and Stromman, 2011; Gold 
and Seuring, 2011). 

The biomass SC configuration has to consider the changes 
of biomass properties along it, to suitably evaluate the biomass 
amounts to be transported and used to produce the final 
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Superscripts and subindices 

or as received basis 

da/ dry and ash free basis 

e electric 

th thermal 

Acronyms 

AWR agricultural woody residues 

BD bulk density 

BM biomass 

CCS carbon capture and storage 

CES centralized energy systems 

d day 

DM dry matter 

EcO economic optimization 

EnvO environmental optimization 

Eff efficiency 

FU functional unit 

FWR forest wood residues 

GHG greenhouse gases 

GIS geographic information system 

IGCC integrated gasification combined cycle 

LCA life cycle assessment 

LCI life cycle inventory 

LCIA life cycle impact assessment 

LHV lower heating value 

MC moisture content 

MILP mixed integer linear program 

MINLP mixed integer non-linear program 
moMILP multi-ojective mixed integer linear program¬ 
ming 

NPV net present value 

OIL bio-oil 

O&M operation and maintenance 

PC pulverized coal 

PEL pellet 

pts points (from LCA methodology, Impact 2002+) 

SC supply chain 

SCN scenario 

SLU bio-oil 

SOTA state-of-the-art 

STN state task network 

TOP pellet of torrefied biomass 

TOR torrefied biomass 

UTM Universal Transverse Mercator 

yr year 

Mathematical formulation: Indices 

a mid point environmental impact categories 

e suppliers 

f,f facility locations 

g end point environmental impact categories 

i tasks 

;' equipment technology 

k intervals for piecewise linear approximation of 

economies of scale 
s materials (states) 

t, / planning periods 


Sets 

A g set of midpoint environmental interventions 
that are combined into endpoint damage fac¬ 
tors g 

E rm set of suppliers e that provide raw materials 
Eprod set of suppliers e that provide production ser- 

Etr set of suppliers e that provide transportation 
services 

FP set of materials s that are final products 
I set of tasks i with variable input 

Ij set of tasks i that can be performed in technol- 

°g y j 

) e technology; that is available at supplier e 

Jf technology; that can be installed at location/ 

jfj technologies that can perform task i 

Jstor technologies to perform storage activities 

Mfet set of market locations 

NTr set of production, or non-transport, tasks 
RM set of materials s that are raw materials 
Sstor set of materials/biomass that if stored change 
their properties 

Sup set of supplier locations 

T s set of tasks producing material s 
T s set of tasks consuming material s 

Tr set of distribution tasks 

Parameters 

A s ft (kg) maximum availability of raw material s in 
period t in location/ 

Adapt] (kg) lower bound for the capacity of co¬ 
combustion in plant/ 

Adaptj (kg) upper bound for the capacity of co¬ 
combustion in plant/ 

Dem s f t (MJ) demand of product s at market/in period t 
Cost]’ 0 ® 1 (€) cost of coal supplied to facility/ 
distance^ (km) distance from location/to location/ 
FCFJjf t (€/h) fixed cost per unit of technology; capacity 
at location/in period t 

FE]™' (h) increment of capacity equal to the upper limit 
of interval fe for technology; in facility/ 
i r (adim.) discount rate 

HVj 0aI (MJ/kg) heating value of coal supplied to facility/ 
M (adim.) a big number 

NormFg (adim.) normalizing factor of damage category 
9 

Pricey (€/MJ) price of product s at market/in period t 
Price]™ 1 (€/MJ) investment required for an increment of 
capacity equal to the upper limit of interval fe 
for technology; in facility/ 

Tortuosity (adim.) tortuosity factor 
Water s (adim.) moisture for material s 
Water™** (adim.) maximum moisture for task i per¬ 
formed in equipment j 

Greek symbols 

a s ij (adim.) mass fraction of task i for production of 
material s in equipment;' 

a s ij (adim.) mass fraction of task i for consumption of 
material s in equipment;' 
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Pjf (adim.) minimum utilization rate of technology j 
capacity that is allowed at location/ 
fag (adim.) g end-point damage characterization factor 
for environmental intervention a 
rjjj (adim.) efficiency of task i performed in equipment; 
r^ lant (adim.) efficiency of combustion in plant/ 

Sjjj/ (h/kg) capacity utilization rate of technology j by 
task i whose origin is location/and destination 
location/ 

p^g .,f (€/kg) unitary transportation costs from location/ 
to location/ during period t 
T ijfot («/kg) unitary cost associated with task i performed 
in equipment j from location / and payable to 
external supplier e during period t 
T sfet (€/kg) unitary cost associated with handling the 
inventory of material s in location/and payable 
to external supplier e during period t 
Xest (€/kg) unitary cost of raw material s offered by exter¬ 
nal supplier e in period t 

fijff'a (points/kg) a environmental category impact CF 
for task i performed using technology j receiv¬ 
ing materials from node / and delivering it at 
node/ 

i//^ aI (points/kg) a environmental category impact CF for 
coal supplied to facility/ 

1 /rj. (points/(kg km)) a environmental category impact 
CF for the transportation of a mass unit of 
material over a length unit 

Binary variables 

Vjp. 1 if technology j is installed at location / in 
period t, 0 otherwise 

SOS2 variables 

gjfitt variable to model the economies of scale for 
technology j in facility / at period t as a piece- 
wise linear function 

Continuous variables 

Adaptf installed capacity for co-combustion in plant/ 
DamCgjt (point) normalized endpoint damage g for loca¬ 
tion/in period t 

DamCg C (point) normalized endpoint damage g along 
the whole SC 

EPurchet (€) economic value of purchases executed in 
period t to supplier e 

ESalest (€) economic value of sales executed in period t 
FAssett (€) investment on fixed assets in period t 
FCostt (€) fixed cost in period t 

Fjft (h) total capacity of technology; during period t at 
location/ 

FE(h) capacity increment of technology; at location/ 
during period t 

HV S (MJ/kg) lower heating value for material s 
IC a ft (point) midpoint a environmental impact associ¬ 
ated to site / which rises from activities in 
period t 

Impactj 002 (point) total environmental impact for site/ 
lmpact 2 °° 2 aXl (point) total environmental impact for the 
whole SC 

NPV (€) economic metric, net present value 


Pijjft (kg) specific activity of task i, by using technology 
j during period t, whose origin is location/and 
destination is location/ 

Pro/tt (€) profit achieved in period t 
PVsjjft (kg) input/output material of material s for activity 
of task i with variable input/output, by using 
technology j during period tin location / (this 
must be a production activity) 

Pro/tt (€) profit achieved in period t 
Purch^l (€) amount of money payable to supplier e in 
period t associated with production activities 
Purch ™ (€) amount of money payable to supplier e in 
period t associated with consumption of raw 
materials 

Pure hjj (€) amount of money payable to supplier e in 
period t associated with consumption of trans¬ 
port services 

Sales s jj' t (€) amount of product s sold from location/in 
market / in period t 

S s jt (kg) amount of stock of material s at location / in 
period t 

TCostp" 1 (€) total cost of replaced coal during period t 


bioproduct required by the user (Gold and Seuring, 2011; 
Rentizelas et al., 2009b). Biomass as a commodity should 
meet the adequate quality parameters and standards, while 
promoting the development of the most biomass productive 
areas. Legislation for this purpose is underway and needs 
decision-making tools to support the correct policies and 
appropriate subsidies (e.g. grants, tax incentives, bonuses, 
feed-in tariffs, emissions trading) to establish new markets 
(Dasappa, 2011; Lamers et al., 2012). 

Co-firing biomass with fossil fuels is regarded as an alter¬ 
native of transition toward a carbon free society (Berndes et al., 
2010; Faaij, 2006; Gomez et al., 2010b). Co-combustion or co¬ 
firing can be defined as the simultaneous combustion of two 
or more fuels in the same combustion plant (Berndes et al., 
2010), using at least one type of fossil fuel (Sami et al., 2001). It 
is focused on existing installations, originally performed to 
use 100% fossil fuels to generate from tens to hundreds of 
MW of power (Rodrigues et al., 2007). Reduction of CO2 (taking 
into account a life cycle assessment - LCA), sulfur and nitro¬ 
gen emissions (regarding biomass composition) are directly 
derived from the coal fraction substitution with biomass. Co¬ 
combustion is the most direct application among biomass and 
renewable sources use (Faaij, 2006). It may cost 2-5 times less 
than any other biomass process to produce electricity (Berndes 
et al., 2010). Coal is a primary resource that is predicted to 
last as a major fuel for power generation for decades, there¬ 
fore, efforts are directed toward increasing the plant efficiency 
while decreasing its environmental impact (Giddey et al., 2012; 
Royo et al., 2012). 

The biomass properties and quantity should mimic and 
match as much as possible the replaced fossil fuel. There¬ 
fore, biomass should prevent high moisture content (MC) or, 
low dry matter (DM), low grindability, low bulk density (BD) 
and low lower heating value (LHV). Biomass heterogeneity and 
inherent characteristics, as well as a seasonal production, may 
be overcome using biomass densification, stabilization and an 
appropriate storage planning (Puigjaner, 2011, Chapter 2). For 
those reasons pre-treated or “enhanced” biomass is crucial for 
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the development of the bioenergy sector. Enhanced biomass, 
as known nowadays, ranges from solids like pellets, torrefied 
biomass or torrefied pellets to liquids like pyrolysis oil or bio¬ 
oil (Uslu et al., 2008). 

Several studies evaluate the effect of the coal replacement 
fraction on power plant efficiency, showing that there is a 
range for coal substitution without either efficiency penalty 
nor technical problems in the plant. The work by Dai et al. 
(2008) exhaustively revise the issues that may arise from 
biomass use with coal, ranging from slagging and fouling, to 
catalyst poisoning during flue gas cleaning. The work by Liu 
and Balasubramanian (2014) propose a pre-treatment step to 
avoid ashes issues. The paper by Gomez et al. (2010b) pro¬ 
poses coal substitution between 5% and 10% on energy basis, 
being the 5% the minimum fraction that may produce a valu¬ 
able change in plant performance. The affordable shares on 
an energy basis according to Berndes et al. (2010) are 15% 
for fluidized bed boilers and 10% for pulverized coal (PC) and 
grate-fired boilers. 

The energy efficiency of a coal fired plant oscillates between 
30% and 45% on a LHV basis, depending on the boiler technol¬ 
ogy, subcritical or supercritical water, and the antiquity of the 
plant (Berndes et al., 2010), being 36% a common average value. 
In Van Den Broek et al. (1996), the efficiency ranges from 39% 
to 44% for co-firing in PC power plants, up from 100 MW. The 
work by Gomez et al. (2010a) uses a value of 38% as co-firing 
efficiency. Some other works, such as Sami et al. (2001), shows 
the synergetic effect between coal and biomass. It means that, 
not only the efficiency of the coal combustion plants can be 
maintained, but it may increase when using biomass. 

This work deals with the strategic and tactical levels of 
the SC problem at the design step. As in a classical SC 
problem, processing sites, distributors, transportation, ware¬ 
houses, raw material suppliers and customer markets are 
identified. Specifically important for bioenergy problems is 
that raw materials can be highly distributed, with the excep¬ 
tion of energy crops, so logistic issues have to be properly 
considered (Caputo et al., 2005); biomass production, stor¬ 
age and transportation may represent a significant fraction 
of the whole SC cost (Panichelli and Gnansounou, 2008) and 
of the environmental impact. The approach used here for¬ 
mulates a large-scale mixed integer linear program (MILP) 
that captures the relevant information for each facility and 
flow (Graves and Tomlin, 2003). The bio-based SC problem for¬ 
mulation and resolution face new challenges to capture (i) 
biomass properties change, (ii) different types of biomass com¬ 
bination, (iii) biomass seasonality, (iv) demand seasonality, (v) 
multi-site production of biomass and the subsequent, and (vi) 
size of the problem. This leads to a highly computationally 
demanding model. The starting point is the mathematical 
formulation developed in Puigjaner (2011, Chapter 3), and 
Perez-Fortes et al. (2012): it has been adapted and extended 
to consider a new real complex scenario, a whole country, 
contemplating a substantial number of sites and variables. 
A superstructure of pre-treatments and monthly periods 
are considered. As a result, the model allow to (i) locate 
and identify the most appropriate pre-treatment technolo¬ 
gies and feedstock suppliers to replace a certain amount of 
coal in the already existing Spanish park of thermal power 
plants, (ii) to evaluate matter flows between nodes and ech¬ 
elons and (iii) to give a preliminary approach to the monthly 
sequence, i.e. inputs and outputs by month, production rates 
and inventory, for all the units and biomass production 
sites. 


2. Literature review and the SC 

The biomass SC problem may be addressed using a wide range 
of decision-maker outlooks. As example, Caputo et al. (2005) 
evaluate the net present value (NPV) of 100% biomass projects, 
focusing on transportation. Bowling and Dale (2011) look for 
an optimal SC for a biorefinery, considering overall sales and 
costs optimization to discern between a distributed or central¬ 
ized structure with special attention on transportation costs. 
Ayoub et al. (2009) focus on costs and environmental impact 
through emissions to air, water pollutants and solid wastes. 
Damen and Faaij (2006) perform a life cycle inventory to com¬ 
pare co-combustion and combustion of only coal and Perry 
and Rosillo-Calle (2008) focus on C0 2 emissions along the 
whole SC. A more recent work from Mele et al. (2011) combines 
the use of mathematical programming with LCA, to perform 
a multi-objective optimization based on the NPV and the LCA, 
to produce bioethanol from sugar cane in Argentina. Environ¬ 
mental evaluations often take into account a LCA (Cherubini 
and Stromman, 2011). Other attempts have been recently done 
to add the social criterion to the economic and environmental 
points of view, as the creation of places of job (Perez-Fortes 
et al., 2012; You et al., 2012). 

The most state-of-the-art works combine multi-objective 
optimization and mathematical programming (MILP, mixed 
integer non-linear program, MINLP, with and without uncer¬ 
tainty and risk consideration) or scenario-based optimization 
with geographic information systems (GIS) for spatial data 
analysis. The literature review from An et al. (2011a) exposes 
that bioenergy is approaching to an important grow and needs 
to integrate strategic, tactical and operational decisions (i.e. 
the operations research point of view) to enhance and secure 
their viability, even if planning models have not been fully 
required (and therefore, developed) yet. The review by Scott 
et al. (2012) points out that the most used method for multi¬ 
criteria decision-making in bioenergy systems is optimization 
using few alternatives (i.e. a limited number of scenarios), 
followed by mathematical optimization (i.e. coverage of all 
possible options). Usually, the performed case studies consider 
one single biomass type. An et al. (2011b) demonstrate the 
economic viability of a biofuel SC (lignocellulosic biomass to 
bioethanol) in Texas through an strategic-tactical mathemati¬ 
cal model. The works by Zamboni et al. (2009,2009) formulate a 
MILP model to minimize operating costs and GHG emissions of 
a biofuel SC. The second work performs a multi-objective opti¬ 
mization by taking into account the whole SC, from biomass 
cultivation to fuel distribution. The region under study is dis¬ 
cretized into a grid of square regions to map the network 
and to estimate the biomass cultivation potential for each 
region. Using as starting point the formulation developed in 
Zamboni et al. (2009), Akgul et al. (2011) also minimize costs, 
while decreasing calculation time including spatial restric¬ 
tions through a “neighborhood" flow limitation. The starting 
point in Giarola et al. (2012) is Zamboni et al. (2009). The 
approach here is extended to consider more than one type of 
raw biomass to produce bioethanol, and deals with biomass 
availability limitation mainly due to food and land contro¬ 
versies. Kin et al. (2011) develop a MILP model to decide the 
location and capacity of fast pyrolysis and Fischer-Tropsch 
plants to produce biofuel from distributed biomass sources. A 
MILP determines optimal sizes and locations of biofuel plants 
that use an energy crop as raw material for already existing 
gas stations in Austria, in Leduc et al. (2008). The optimiza¬ 
tion deals with investment and operating costs. The paper by 
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Fig. 1 - A generic bio-based SC that produces and distributes electricity. Transportation and storage may occur when 


Elia et al. (2011) uses a MILP multi-objective optimization to 
propose a sustainable SC for a novel hybrid concept for a plant 
that produces liquid fuels, i.e. gasoline, diesel and kerosene, 
using coal, biomass or natural gas as main feedstocks. The 
authors conclude that the hybrid configuration (fossil fuels 
- biomass) allows for competitive biofuels prices. The work 
by Tavares et al. (2011) use GIS to identify the most suitable 
locations for municipal solid waste plants, according to socio¬ 
economic, technical and environmental criteria. Panichelli 
and Gnansounou (2008) solve the location-allocation problem 
of two torrefaction plants and a gasification unit, by means of 
a costs criterion, for a specific region in the north of Spain. 
The exact spatial distribution of biomass waste resources 
roads map are used. Zhang et al. (2011) adopt a GIS-based 
method to find the best location for a biofuel plant, based 
on the distributed nature of the woody biomass sources and 
the associated transportation costs. The methodology con¬ 
sists of two steps: location of site candidates and selection of 
the cost-optimal area. The paper by Chiueh et al. (2012) uses 
GIS to identify the biomass waste resources and the specific 
transportation routes for a co-combustion retrofitting prob¬ 
lem, evaluating different levels of pre-treatment (torrefaction) 
centralization. 

Different works can be found that go in depth into a specific 
SC echelon. Fig. 1 represents a generic bio-based SC to produce 
and distribute electricity, with multiple biomass providers and 
multiple biomass processing plants: (i) feedstock production 
(growing, harvesting and collection), or waste generation, (ii) 
biomass pre-treatment, (iii) storage, (iv) biomass treatment, 
(v) electricity distribution and (vi) electricity consumption. 
Storage and transportation can be needed in any moment 
along the SC, depending on the specific distances to be 
covered, the supply and the demand characteristics. Never¬ 
theless, it is more convenient to store and transport biomass 
in an upgraded state, i.e. after a pre-treatment, to avoid costs 
and non-desired effects. Transportation and storage are the 
most critical echelons for a bio-based SC (Caputo et al., 2005; 
Iakovou et al., 2010; Rentizelas et al., 2009b). The electric¬ 
ity distribution level may include the desired level of detail 
of a traditional grid, a smart grid or a microgrid, depending 
on each case study. Storage is imperative when considering 
seasonal biomass supplies and electricity demand. Biomass 
different origins, high MC, low BD, low LHV, and biomass 
degradation during storage, lead to a needed improvement 
and homogenizing step to optimize transportation, handling 
and processing, to obtain a “standard” product. Therefore, a 
bio-based SC may consider changes on biomass characteris¬ 
tics into each SC echelon. 

Covering all biomass states and pre-treatments, the means 
of transport used can be truck, trailer, train, ship or pipeline. 


The study by Cherubini and Stromman (2011) reviews the 
works performed in the field of LCA for bio-based systems, 
showing a general reduction of GHG emissions and of fos¬ 
sil fuel consumption when using biomass instead of fossil 
fuels. Even if the LCA follows a well-established methodology, 
the selection of the functional unit (FU) changes according to 
authors criteria (for instance, power produced or energy/mass 
flowrates of biomass introduced into the system) and it makes 
difficult the comparison among different works. Furthermore 
many LCA’s have been conceived under different assumptions, 
hence their comparison is also difficult. 

According to the reviewed works, no other work has taken 
into account a superstructure of pre-treatments solved by a 
MILP formulation that considers simultaneously the change 
of biomass properties along the different echelons, to evalu¬ 
ate large scale retrofitting co-combustion plants projects. The 
gaps to cover to successfully implement bionenergy projects 
are: (i) integration of decision levels, (ii) multi-objective opti¬ 
mization, (iii) biomass properties (mainly MC and BD) change 
along the SC, (iv) pre-treatments evaluation, (v) seasonal 
raw matter and demand, and (vii) limitation of calculation 
time. According to that, the approach of this work pro¬ 
poses a sustainable SC by supporting the decision-making 
task regarding production/distribution and storage network, 
selection of most appropriate technologies, units sizing, allo¬ 
cation of matter flows, allocation of products and monthly 
planning of activities, while providing costs and environ¬ 
mental impact shares. The co-combustion echelon in the 
mathematical formulation is modeled as a transformation 
coefficient to estimate electricity production, i.e. the pro¬ 
cess efficiency on LHV ar basis, 36% as a conservative value 
(Berndes et al., 2010). The problem assumes a maximum 
coal substitution of 15% in thermal basis and by power 
plant. It deals with monthly periods for a given amount 
of identical years. It considers solid and liquid biomass 
transportation. The developed tool is used for a specific 
case study, but can be adapted to handle any other type 
of case, with wider (international) or narrower (regional) 
frontiers. 

The article is organized as follows. The superstructure of 
biomass pre-treatment technologies is characterized, together 
with biomass storage, as echelons that change biomass prop¬ 
erties. After this section, the major features of the approach 
are described. Then, the problem statement sums up the main 
input and output data. The mathematical model part depicts 
the current model from the previous formulation developed 
(Perez-Fortes et al., 2012). Following that, the paper focuses 
on the specific case study, and finally on the results and dis¬ 
cussion of the economic and environmental results of the 
scenarios performed. 
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3. Biomass storage and superstructure of 
biomass pre-treatments 

The presence of pre-treatment technologies is needed to 
enhance biomass properties to convert biomass into a com¬ 
modity in energy projects. Thus, pre-treatment requirement 
is one of the drawbacks of using biomass. The work by Bals 
and Dale (2012) remarks that as the pre-treatment operations 
increase considerably the cost of the final biomass product, 
it is crucial to determine if the tasks of homogenization and 
densification are needed to have the desired results. As alter¬ 
native, it can be determined if the project can be locally 
performed, that is the promotion of distributed projects rather 
than CES. Nowadays, biorefinery installations have not been 
widely established since there is a lack of sustainable SC 
for lignocellulosic feedstock. Biorefineries (centralized) need 
distributed pre-treatment installations. However, current pre¬ 
treatment technologies are expensive and inefficient, reasons 
why this field needs more research and development (Kurian 
et al., 2013). Wu et al. (2010) put into manifest that a bioslurry- 
based SC would be competitive if more than 1500 dry tons per 
day are being transported. On the other side, for less than 500 
dry tons processed per day, a conventional SC is favored. 

Fig. 1 depicts five main steps, that modify biomass 
characteristics: biomass collection, pre-treatment, storage, 
transportation and treatment. Pre-treatment and storage ech¬ 
elons may change biomass characteristics, as described in the 
following lines. Torrefaction and fast pyrolysis are the state- 
of-the-art technologies considered, concretely fast pyrolysis 
in an early stage of development, while the well-known pel¬ 
letization is also included by itself, or in combination with 
torrefaction, as pre-treatment options. 

• Storage. This echelon may change MC, DM and LHV due 
to microbiological degradation. It is mainly affected by the 
temperature of the stored biomass, biomass properties and 
state (solid or liquid). The type of storage depends on the 
class of biomass to be stored, i.e. (i) open-air storage in trop¬ 
ical summer or during hot periods for raw solid biomass, 
which even contributes to biomass drying, (ii) silos and 
bunkers (closed storage) for pre-treated biomass. Fermenta¬ 
tion can develop hot spots, which can set on fire the biomass 
piles, so storage facilities must be carefully designed. Liq¬ 
uid biomass, like bio-oil, rarely changes its properties while 
stored if convenient conditions are set. To alleviate degra¬ 
dation, the stored biomass should be as homogeneous as 
possible and have a MC usually under 20% on a mass basis 
(Rentizelas et al., 2009a, b). 

• Torrefaction. It is a thermal process that takes place at 
relatively low temperatures, 225-300 °C, and uses a low 
heating rate 50°C/min, performed at atmospheric pres¬ 
sure in an inert atmosphere (Prins et al., 2006a). The solid 
product is a uniform solid with lower MC and higher LHV 
than the raw material, as well as highly hydrophobic Uslu 
et al. (2008). The by-products obtained are a condensable 
liquid and a non-condensable gas. It is a relatively new 
technique applied to biomass, but the benefits have been 
already perceived via higher BD values and better grind- 
ability (Prins et al., 2006b; Uslu et al., 2008) that enhance 
biomass gasification efficiency (Couhert et al., 2009; Prins 
et al., 2006a), produce a less reactive char with steam after 
gasification (Couhert et al., 2009) and limit biological degra¬ 
dation (Bergman et al., 2005). The paper by Deng et al. (2009) 
evaluates the performance of torrefied agricultural residues 


in co-gasification with coal: torrefied biomass is more simi¬ 
lar to coal than raw biomass, as well as more homogeneous. 

• Pelletization. It is a physical process that dries, mills, densi¬ 
ties and homogenizes biomass, frequently using only lignin 
as binding agent. The BD is higher, biological degrada¬ 
tion is restricted, and the matter is easier to handle than 
raw biomass, with the subsequent benefits in transporta¬ 
tion and storage. Nevertheless, pellets, which are cylinders 
of 6-8 mm diameter, may acquire moisture, have a lower 
mechanical resistance toward crushing and may form dust. 
The market share of pellets is growing and this is the most 
common form of biomass commercialized. Wood is the pri¬ 
mary raw material used (Sultana et al., 2010; Maciejewska 
et al., 2006). 

• Pelletization of torrefied biomass (TOP pellets). The combi¬ 
nation of both techniques can overcome the drawbacks of 
torrefied biomass and pellets separately, with subsequent 
advantages in transportation and storage. TOP pellets have 
low MC with limited predisposition to uptake humidity, 
biological degradation is almost completely inhibited, den¬ 
sity and calorific value are higher. Moreover, TOP pellets 
allow a wider range of raw biomass types and yield simi¬ 
lar physical properties to the final product. Pelletization of 
torrefied biomass is less energy consuming than pelletiza¬ 
tion of untreated biomass (Uslu et al., 2008; Maciejewska 
et al., 2006). According to Uslu et al. (2008), TOP pellets are 
the most suitable option for international trade. 

• Fast pyrolysis. Pyrolysis is a thermal process that decom¬ 
poses biomass at temperatures in the range of 400-800 °C in 
an inert atmosphere. Slow and fast pyrolysis differs in the 
heating rate, and therefore, in the proportion of char, gas 
and liquid products obtained. Products range also depends 
on the biomass type and reaction parameters (Uslu et al., 
2008). Fast pyrolysis takes place at 450-550°C and counts 
with heating rates as high as 60,000 °C/min. Its main prod¬ 
uct is in liquid phase, the so-called bio-oil or pyrolysis oil. 
It is a mixture of 70% oxygenated organics and 30% water, 
on a mass basis (Maciejewska et al., 2006). Bioslurry is the 
mixture of bio-oil and biochar, which is the solid product 
from the pyrolysis process, milled into fine particles. Hence, 
the efficiency of the process is higher. As a liquid fuel, trans¬ 
portation and handling can be easier and cheaper (Abdullah 
et al., 2010; Wu et al., 2010). 

The different biomass properties transformation 
coefficients for each technology and echelon are sum¬ 
marized in Section 7, in accordance with linear models for 
MC, DM, BD and/or LHV modification. Table 5, in Appendix 
A, compiles the biomass properties that change into each SC 
echelon considered. 

4. Problem statement 

The bio-based SC configuration and operation should be tailor- 
made. The approach described in this work can be adapted to 
any other decision-maker needs, for example, the fraction of 
coal to be replaced, different technologies share and biomass 
types, with the subsequent differences on seasonality and 
intrinsic properties. 

The bio-based SC assumes fixed located sites for biomass 
collection, as well as with established points for demand, elec¬ 
tricity generation plants that request a specific amount of inlet 
energy supplied by biomass waste. A list of potential locations 
for pre-treatment sites and distribution centers or storage is 
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proposed taking into account the highest biomass waste pro¬ 
duction sites. Processing sites are limited by minimum and 
maximum capacities. They can be supplied by more than 
one raw matter producer. Investment costs take into account 
economies of scale. Distribution centers and transportation 
activities are modeled by considering upper and lower bounds 
based on their biomass handling capacity. Distribution centers 
can be supplied from more than one biomass producing site 
or pre-treatment plant. A material flow between facilities may 
appear if selecting such flow improves the performance of the 
SC. The market demand for energy can be satisfied by more 
than one site. 

Inputs: 

1. Process data 

• Amount of raw biomass waste that is available. Matter 
characteristics and seasonality. 

• SC echelons and inlet biomass properties required, i.e. 
conditions of humidity, shape and/or energy content, for 
the specified activities. A set of pre-treatment technolo¬ 
gies, types of storage and means of transportation. 

• Activities (pre-treatment, storage, transportation and 
treatment) efficiencies, to linearly model each process 
and determine the outlet matter characteristics in terms 
of MC, DM, BD and LHV. Utilities consumption and life¬ 
time for each one of the activities. 

• A set of matter states that quantifies MC, DM, BD and 
LHV for each flow of mass between activities. 

• A set of demands, i.e. amounts of seasonal energy 
required by the co-combustion plants. 

• A set of providers, intermediates and consumers loca¬ 
tions. 

• Time period, planning horizon, project lifetime and aver¬ 
age working hours for a specific interval of time. 

2. Economic data 

• Investment, fixed and variable costs associated to all the 
technological options involved: pre-treatments, storage 
and co-combustion plants adaptation for biomass firing. 
Unit transportation costs per km and volume of biomass 
to move. 

• Base escale and the associated economies of scale for 
technologies capacity. 

• Price of raw biomass, coal, utilities, consumables and 
electricity. 

• Interest rate. 

3. Environmental data 

• Environmental interventions of the biomass waste and 
coal production. 

• Environmental interventions of the electricity produc¬ 
tion from coal. 

• Activities environmental interventions. 

• Utilities environmental interventions. 

Outputs: 

1. The network structure that shows the optimal level of 
centralization/de-centralization: selected pre-processing 
units with their capacities and locations; types of trans¬ 
portation, storage sites and dimensions, percentage of 
substitution of coal for each power plant. Connections 
among sites. 

2. Feedstock utilization and schedule, i.e. suppliers operation 
per month. 


3. Flows of biomass among sites at each time period. 

4. Inventory levels at each time period. 

5. Breakdown of costs, revenues, cost saving and investment, 
and environmental impacts per echelon or equipment. 
Impact categories share per each one of the options con¬ 
sidered, biomass usage vs. coal replaced. Economic and 
environmental criteria values for the whole time horizon. 

The multi-objective optimization integrates economic and 
environmental criteria, taking into account the NPV, in €, 
and a LCA in Impact 2002+points (pts). The assumed objec¬ 
tive functions quantify the difference between the only-coal 
with the coal-biomass SC. 

• Economic optimization. The NPV evaluates the investment 
and the operational costs belonging to the biomass SC and 
the feedstock savings related to the SC with only coal (dif¬ 
ference between the cost of the coal replaced and the waste 
used), at the same level of incomes. 

• Environmental optimization. The LCA expresses the envi¬ 
ronmental impact difference between 100% coal power 
plants and plants that replace a fraction of feedstock by 
biomass. 

5. Mathematical model 

The mathematical model of the biomass SC considers dif¬ 
ferent decisions among them: (i) activation of SC nodes and 
links among them (a node can be represented by processing 
or distribution activities), (ii) the capacity to carry out each 
period, (iii) their capacity utilization levels, (iv) selection of the 
most appropriate pre-treatment technologies, (v) transporta¬ 
tion links, (vi) fraction of coal replaced by biomass, and (vii) 
the amount of matter moved along the different SC nodes. As 
performance indicators, the overall NPV as well as the envi¬ 
ronmental impact savings associated with the whole SC are 
utilized. 

The resulting model is solved using a multi-objective MILP 
approach, which allows to assess the trade-off between the 
environmental damage categories and the economic indica¬ 
tor. The variables and constraints of the model are classified 
into three groups: (i) process operations constraints given by 
the SC topology (the so-called design-planning model), (ii) the 
economic metric and (iii) the environmental metric formula¬ 
tion. They are described in the following paragraphs. 

Fig. 1 is a general outline of the different components of 
the bio-based SC model: sourcing, pre-treatment, final prod¬ 
uct processing and distribution. The SC is defined as a number 
of potential locations where processing sites or distribution 
centers, or both of them can be installed. Suppliers are at 
fixed locations where biomass is available. The product can be 
processed at several sites. The properties of raw biomass may 
be modified/improved by means of the pre-treatment units so 
as to allow intermediates to meet the characteristics required 
by subsequent steps in the SC. Pre-treatments may be also 
convenient to induce savings in transportation costs. 

5.1. Design-planning model 

The design-planning model selected to deal with the biomass 
SC network is adapted from the works of Lainez-Aguirre et al. 
(2009), Perez-Fortes et al. (2012). This model translates the 
State-Task-Network (STN) formulation (Kondili et al., 1993), 
which is a widely known approach for scheduling, to the 



1546 


CHEMICAL ENGINEERING RESEARCH . 


DESIGN 92 ( 2014 ) 1539-1562 


SC context. The main advantage of this formulation is that 
it collects all the SC nodes activities information through a 
single variable set. This feature eases the economic and envi¬ 
ronmental metrics formulation as well as it facilitates the 
consideration of the different pre-treatment and their possible 
combination to obtain the required outputs. The SC material 
balances can be modeled by means of a single equations set 
for all materials and facilities. This is possible since processing 
nodes and distribution centers, as well as final products, raw 
material and intermediates are handled indistinctively. The 
most relevant variable of the model is P ; p, t [kg], which rep¬ 
resents the magnitude of particular task i, performed using 
technology j during period t, whose origin is location / and 
destination is location /. In the case of production activities, 
they must receive and deliver material within the same site 
(Pjjjjt). In contrast, in a distribution activity, facilities / and / 
must be different. This mathematical formulation assumes 
that an activity consumes and produces certain materials 
with determined properties and can be performed in differ¬ 
ent equipments. By using the activities as the core of the 
formulation rather than products-materials, it results in a 
flexible formulation that can easily incorporate new technolo¬ 
gies/processes. The equations comprising this formulation are 
described in the following paragraphs. 

The mass balance must be satisfied at each node of the 
network. The expression for the mass balance for each type of 
material s (i.e. raw material, pre-treated biomass, or final prod¬ 
uct) processed at each potential site/in every time period t is 
presented in Eq. (1). Parameter a s y [dimensionless] is defined 
as the mass fraction of material s that is produced by task i 
using technology;. T s is a set that refers to tasks that produce 
s, while a s jj and T s are associated with tasks that consume s. 

S sft~ S sft-l = Y 

S' 

-EE E ^p im vsj, t (i) 

S' ieTsMiH?/) 

The change in inventory of material s for consecutive plan¬ 
ning periods and for each location is given by the difference 
between the amount of material produced by those tasks 
belonging to set T s and the amount consumed by those tasks 
included in set T s . The model assumes that process parame¬ 
ters such as conversions, separation factors or temperatures, 
are fixed for each activity in order to enforce the linearity of 
the problem. In this sense, parameters a sij [dimensionless] and 
a S y [dimensionless] provide the “recipe” for a specific activity. 
Nevertheless, there are activities (I) for which it is desirable to 
let the model specify the mixture of inputs in order to achieve 
a given value of a specific biomass property, for instance, a 
specific moisture content or calorific value. For such activities, 
the combination of feedstocks and, therefore, the proportion 
of each feedstock is variable. In order to take into account such 
activities, the mass balance is modified as shown in Eq. (2). 
Note that Eq. (1) is a particular case of Eq. (2). 

S sJ t - Ss/t-i = EE E asi J P <ifft ~ EE E 

S' ieTsjeftnJf,) S' ieT s ;e(JiO?/) 

+ E E Pu ^ - E E Pv ^ Vs ’ f ’ t 

ie(T s nl)je(Jjnj/) ie(T s ni)MiOj/) 

(2) 


With regard to the variables Pyjft [kg] and Pu si p [kg], note 
that the former is used in the mass balance (Eq. (2)) coupled 
with the parameter a sjj - or a s y which specify a fixed proportion 
of material produced or consumed for a task i. On the other 
hand, Pu s ;p is modeling flexible tasks which allow the propor¬ 
tion of the material produced or consumed to vary so as to 
provide alternatives for biomass mixing to the model. For this 
reason, Pu s jp is not multiplied by such parameters in Eq. (2). 

For the activities for which biomass properties are fixed, 
the energy balance is satisfied directly by the definition of the 
streams. However, it is necessary to verify that the energy 
balance is satisfied for the flexible activities (I). The energy 
balance for such activities is represented by Eq. (3). Here, HV S 
[MJ/kg] is the heating value of material sand r/y is the efficiency 
of process i performed in equipment j. Each different type of 
biomass has a different heating value. Therefore, a specific 
activity changes the heating value of the output stream if (i) 
it is a pre-treatment task that modifies explicitly the calorific 
value of the biomass, or (ii) it is a task whose main objective is 
the change of shape, but it is receiving a mixture of biomasses 

mj^HV s Pv si j ft = ^HV s Pv sljft Vi e I ,j,f t (3) 

^T s se f s 

Another important extension of this model is the consid¬ 
eration that storage is capable of changing biomass properties 
(e.g. MC, DM, LHV). In order to do so, storage must be consid¬ 
ered as an actual activity. Subsets J s tor and S st0 r will represent 
the storage “equipment” and those materials that when kept 
in storage change their properties, respectively. Notice that the 
mass balance has been decomposed in this case so as to deal 
with the one period delay necessary for the properties change 
to occur. As expressed by Eqs. (4) and (5), a storage activity 
places inventory in the current period t and takes inventory 
from the previous period t — 1. 

S *=EE E «■> vse w.t < 4 ) 

S’ 'eTsMin/f'OJstod 

S *- 1= EE E Vs e Sstor,/, t (5) 

S' fefsMrf/nW) 

In order to keep the model linear Eqs. (2) and (3) are 
limited to combine different materials in proportions which 
are defined by the optimization just for the activities right 
before the energy generation. The states feeding any other 
activity and their corresponding proportions and properties 
(e.g. MC, DM, BD, LHV) must be defined a priori. 

Eqs. (6) and (7) represent the temporal change in the equip¬ 
ment technology installed in a potential facility location. We 
will consider economies of scale by using a piecewise linear 
approximation in K different intervals and a so-called SOS2 
variable type (fp,). Such variables are positive and at most 
two consecutive variables are non-zero. FEj™ 1 [h] is the limit 
of capacity for interval k. Vp is a binary variable indicating 
whether or not capacity of technology j is expanded at site / 
in period t. This formulation will be recalled in the economic 
performance metric section for computing the investments 
associated with capacity expansions. Eq. (8) is used for the 
total capacity Fp [h] bookkeeping taking into account the 
capacity expansion during the planning period t (FEp [h]) for 
equipment technology j in facility/. This equation considers 
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the case of the initial design of a SC (FE^ 0 = 0) as well as a SC 


retrofit scenario (FEp =f= 0). 


E^ tFE ^ it = F % 

(6) 

ZM = % Vf,jejf.t 

(7) 

F jft = Fjft^ + FEj ft Vf, jejf.t 

(8) 


Eq. (9) is used to ensure an utilization greater than or equal 
to a minimum value established by the decision-maker and 
that the utilized capacity is lower than or equal to the avail¬ 
able one. Parameter ftjf [dimensionless] defines a minimum 
utilization rate of technology j in site / as a proportion of the 
total available capacity. Parameter 9^ [h/kg] represents the 
resource utilization factor. This is the capacity utilization rate, 
in terms of capacity units (e.g. machine-hours), of technology 
j by task i whose origin is location / and destination location 

f- 

s EEVVt < %-1 Vf, jejf.t (9) 

/' ielj 

The capacity is expressed as equipment j available time 
during one planning period, then 9^ represents the time 
required to perform task i in equipment j per unit of pro¬ 
duced material. Thus, once operation times are determined, 
this parameter can be readily approximated. 

Eq. (10) guarantees that the amount of raw biomass s pur¬ 
chased from site / at each time period t is lower than an upper 
bound given by suppliers availability A s ^ [kg] (e.g. seasonality, 
residues availability in a specific region). 

^ Vs e RM ’f e Su P > 1 ( io ) 

/' ieT s Mi 

Adaptf is a variable which represents the capacity of a co¬ 
combustion plant to process biomass fuels, while r^ ,ant is the 
efficiency of plant/. Eq. (11) expresses the feasible limits for the 
modification of a combustion plant to process biomass fuels. 
Eq. (13) expresses that inlet energy flows must be lower than 
the available capacity at the co-combustion plant/to process 
biomass fuels. 

By Eq. (12) sales of final product s e FP carried out from facil¬ 
ity location/ to market/e < Mfet are estimated. Eq. (14) is used 
to express that the demand can be partially satisfied due to 
biomass production or supplier capacity limitations: the sales 
of product s carried out in market / during the time period t 
could be less than or equal to the demand. 

Adaptf < Adaptf < Adapty V/ e Mfet (11) 


Sa,es sfft = 9f ant HV s ( £ X)P iif J 

V e (T= nTr )je(JiO?f) / 

xVseFPJe Mfet, /' i Mfet, t (12) 


^2 ^2 S Q l es s /'jt - Adaptf Vf e Mfet, t 

(13) 

seFP/VMfet 


yy Saies..f/ft < Dem s fl Vs e FP,/e Mfet, t 

(14) 


ft Mfet 

It is important to emphasize that for this sort of networks 
the final product is the energy delivered to the different 
regions (Mfet). Thus, for those tasks that carry out the energy 
generation their corresponding parameter a si j is a function of 
the heating value of the input materials and the efficiency of 
the equipment. 

5.2. Economic formulation 

The expressions required to compute the operating revenue, 
the operation costs, the total capital investment, and NPV are 
next described. 

The operating revenue is expressed in Eq. (15) as the prod¬ 
uct sales during period t. 

ESalest = ^2 Sales s fff t Price s f t Vt (15) 

seFP/eMfetfy(MfetuSup) 

The operating costs include fixed and variable costs: Eq. 
(16) describes the total fixed costs of operating the SC network. 
FCFJjfl [€/h] is the fixed unit capacity cost of using technology 
j at site/. 

FCostt = J2 J2 FCF} ifi F >fi Vt (16) 

/ ^(MfetuSup) jejy 

In turn, as variable costs, the cost of purchases from sup¬ 
plier e, includes raw material procurement, transportation and 
production resources, as shown in Eq. (17). The purchases of 
raw materials (Purch™ [€]) made to supplier e are evaluated 
in Eq. (18). We will assume a different supplier for each com¬ 
ponent of the variable costs. This assumption can be easily 
relaxed to account for the specific characteristics of the prob¬ 
lem being dealt with. The variable xest [€/kg] represents the 
cost associated with raw material s purchased to supplier e. 
Transportation and production variable costs are determined 
by Eqs. (19) and (20), respectively. p^ t [€/kg] denotes the e 
provider unit transportation cost associated with material dis¬ 
tribution from location / to location / during period t. 
[€/kg] represents the unit production cost associated with 
performing task i using technology], whereas [€/kg] rep¬ 
resents the unit inventory costs of material s storage at site 
/. The parameter entails similar assumptions to the ones 
considered for a s ,j and a S y, namely, the amount of utilities and 


labor required by an activity are proportional to the amount 
of material processed. 

EPurch e t = Purch™ + Purchf t + Purch^ od Ve, t 

(17) 

Purch™ = ve e E™, t 

seRMfeF e ie f s jeJi 

(18) 

Purchft = ^ 'E'E P wM}f't VeeEtr.t 

* eTr jtli n Je f ? 

(19) 
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?urc\C d = EE E p yfft^fet 

f ifTrjeVirff) 

+ E E Ve6 W ( 20 ) 

s /^(SupUMkt) 

The total capital investment on fixed assets is calculated 
by means of Eq. (21). These equations include the investment 
made to expand the technology’s capacity j in facility site/in 
period t. Recall that economies of scale for technologies capac¬ 
ity is considered using a piecewise linear approximation in K 
intervals. Here, Price]™ 1 [€] is the investment for a capacity 
expansion equal to the limit of interval fe (FEj™ 1 [h]). 

FAssefi = EEE Price ^V vt ( 21 ) 

/ j k 

Eq. (22) calculates the profit in period t, as operating rev¬ 
enues minus fixed and variable operating costs. We take into 
account the ‘cost’ related to the coal that is replaced by 
biomass which is expressed in Eq. (23). The NPV is calculated 
as in Eq. (24). Note that the NPV only for the biomass- 
based-energy generation and considers the ‘opportunity’ costs 
associated with the fraction of replaced coal. 


Consequently, the value of the mid-point environmental 
impact i'jjff'a associated with transport, is calculated as in 
Eq. (26), where [pts/(kgkm)] represents the a charac¬ 
terization factor of the environmental category impact for 
the transportation of a mass unit of material over a unit of 
length. Note that the impact is assigned to the origin node. 
The environmental impacts associated with production (Eq. 
(25)) or transportation (Eq. (26)), can be calculated by setting 
the indices summation over the corresponding tasks (i e TV or 
i e NTV). A tortuosity factor (Tortuosity [dimensionless]) may be 
employed to correct the estimated (linear) distance between 
nodes. 

= ifJ- a distanceg’Tortuosity Vi s Tr, j e/, a,//' (26) 

Eq. (27) introduces DamC^ [pts] which is a weighted sum 
of all mid-point environmental interventions. They are com¬ 
bined using g end-point damage factors f„g [dimensionless], 
normalized with NormF g factors. Moreover, Eq. (28) computes 
the g normalized end-point damage along the SC (DamCg C 
[pts]). 

DamCgft = NormFg hglCgft Vg,/,t (27) 

ae. A g 


Profit t = ESalest - ^FCosfi + ^^EPurch et ^ Vt 

(22) 

rnn, V-V^V^ C0St f 0 a ' SaleS sf'fi 

feMktf?Mkt s 'J f 

(23) 

NPV V ( PWfitt ~ FASSStt + TCost t C ° al ^ 

vl (l + ir) f J 

(24) 


5.3. Environmental formulation 

The application of the LCA methodology to the SC model 
requires four steps according to ISO documents (ISO14040): 
goal definition and scope, life cycle inventory (LCI), life cycle 
impact assessment (LCIA) and results interpretation. LCI val¬ 
ues are retrieved from Ecoinvent (2006) database using PRe 
Consultants (2004). 

Environmental interventions are translated into metrics 
related to environmental impact as end-points or mid-points 
metrics using characterization factors. Eq. (25) calculates I C a ft 
[pts] which represents the mid-point environmental impact a 
associated with site/, as a consequence of carrying out activi¬ 
ties in period t. In turn, 'l'ijjf a [pts/kg] is the a characterization 
factor of the environmental category impact for task i per¬ 
formed using technology j, receiving materials from node / 
and delivering them to node/. 

%=EEE%«%'t va -f’ 1 < 25 ) 

je/f ieI i S' 


Likewise to a s y and a s jj, the value of i'gg'a is fixed and con¬ 
stant, provided that all environmental impacts are considered 
linearly proportional to the activity performed in the node 
(Pjjjjt [kg]) (Heijungs and Suh, 2002). Environmental impacts 
associated with transportation, have as FU the amount of 
kg of material transported over a given distance (kgkm). 


DamCg C = y^y^ DamC g ft Vg (28) 

f t 

Similarly to the economic formulation, we consider the 
impact that is avoided by replacing the coal by biomass in 
the combustion plants. Eq. (29) expresses the impact of the 
replaced coal. 


Impact Coal = EEEEE 

g aeAgfeMktftfMkt t 


ifrCoai N ormFg£ ag Sales s fft 

r/fLHVf™ 1 


(29) 


The impact of coal considers (i) coal supply, where 
imported coal, moreover, considers the impact of the coverage 
of longer distances, and (ii) electricity generation as a whole, 
including (iii) transportation impacts. Eq. (30) collects the end¬ 
point environmental damages for the whole SC considering 
the savings in terms of impact due to the replacement of coal. 
Biomass SC impact distinguishes between solid and liquid 
transportation. The impact associated with each processing 
unit is proportional to the utilities consumption, i.e. tons of 
diesel consumed or GJ of electricity used. 

Impact™, 2 ^! = Impact Coal - (30) 

/ 3 1 

For further details regarding the operational and environ¬ 
mental formulation, and the specific points distributed by 
activities, the interested reader is referred to Bojarski et al. 
(2009) and Perez-Fortes (2011). 

The overall optimization problem can be posed mathemat¬ 
ically as follows: 

Max | NPV, Impact™ 2 a „ ] 
subject to 
Eqs.(l)-(30); 

Xe {0, l)Jsf+ 
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Fig. 2 - Specific mathematical formulation nomenclature. The scheme points out how to deal with biomass mixtures and 
properties changes in the linear program: combinations of states and tasks. 


where, X denotes the binary variables set, while y corresponds 
to the continuous variable set. Note that NPV and Impact™ 2 aII 
represent the differences between the SC with coal and the 
biomass SC, being negative and positive, respectively. 

6. Major features in the bio-based SC 
approach 

The model intends to elucidate strategic decisions like the 
final distribution network and selection of pre-treatment tech¬ 
nologies, as well as tactical determinations like the planning of 
processes, flows between nodes and allocation of products. It 
contributes by integrating a superstructure of pre-treatments 
and considering biomass properties MC, DM, LHV and BD 
along the SC. The model is flexible to incorporate/update infor¬ 
mation or change the case of study. 

The linearity of the mathematical formulation does not 
allow the multiplication of variables, like the flow of differ¬ 
ent inlet biomass that are combined conveniently to create a 
“new” biomass with different properties MC, DM, LHV and BD. 
This happens in the SC proposed when considering (i) biomass 
storage, (ii) biomass combinations before pre-treatment, and 
(iii) biomass combinations before treatment. To deal with that, 
the current model proposes the definition of specific states 
(s) and tasks (i) that assumes and defines a priori the pos¬ 
sible results of the storage activity or of each one of the 
pre-treatments. See in Fig. 2 an outline of the different ech¬ 
elons and the related mathematical nomenclature. Since the 
co-combustion plants needs a specific amount of energy (flow 
inlet condition), the only important characteristic here is the 
LHV multiplied by the mass flowrate. Moreover, the output of 
the co-combustion plants is not a biomassic state defined by 
its properties. For that reason, this echelon is modeled using 
the so-called flexible units, explained in Section 5. 

This problem assumed a maximum of two months at open 
air storage for raw material, and its subsequent biomass from 
different periods mixtures. No combination between forest 


wood residues (FWR) and agricultural wood residues (AWR) 
is allowed. It creates a total of 55 states along the whole SC. 
Nevertheless, after preliminary model runs and due to the spe¬ 
cific conditions of the problem, the model decides to store only 
pre-treated matter, since its storage does not change biomass 
properties. So, for the reported case studies results, we sim¬ 
plified the number of states to 15. 

The model is solved along a period of 10 years, assuming 
yearly strategical decisions and monthly planning solutions. 
The MILP is implemented in GAMS, and counts with 59,057 
equations, 5,925,346 continuous variables and 589 discrete 
variables. It takes a maximum of 13,971 CPU seconds to reach 
a solution with a integrality gap of 2.5% on a desktop computer 
dual four core Intel Core i7 920, with 2.66 GHz CPU and 24 GB 
RAM with 8 processors, using the MIP solver CPLEX 9.0. 

7. Case study: a BSC located in Spain using 
co-combustion 

The case study is a retrofitting proposal for coal combustion 
power plants in Spain that contemplates the use of biomass to 
replace a fraction of coal. Given a set of biomass collection sites 
and the list of power plants in Spain, the SC model will help 
on the decision-making task by providing solutions for the 
location-allocation problem and the flows of matter between 
sites, while quantifying them in monetary and environmen¬ 
tal terms. Universal Transverse Mercator (UTM) geographic 
coordinate system is used; linear distances between sites are 
calculated and corrected by a tortuosity factor. The bound¬ 
aries of the problem are from cradle-to-gate: the final echelon 
is the co-combustion plant. The distribution and use of elec¬ 
tricity will be the same irrespective on their source, thus 
selecting a boundary limit at the co-combustion plant seems 
reasonable. On the other hand the environmental interven¬ 
tions associated with the co-firing are not distinguished and 
the environmental impact of co-firing is mainly associated 
with the use of coal. Several issues might arise when co-firing, 
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Table 1 - Feedstock characteristics (Gomez et al., 2010a). 


Biomass MC (wt%) LHV ar BD (kg/m 3 ) Yearly Adjusted Seasonality Cost (€/ton) 

waste (MJ/kg) available (kton) availability (kton) 

FWR 30 12.5 140 7748 1162 None 56 

AWR 40 10.8 100 3883 2718 Summer and winter 52 


see for example (Dai et al., 2008), one of them is ashes forma¬ 
tion, which is increased by biomass use. In this approach we 
assume the use of biomass while co-firing does not changes 
dramatically the environmental interventions of coal power 
production. In this sense these interventions are not signifi¬ 
cant when plant operates at low co-firing ratios, and/or when 
high-quality biomass is used. The currency used is Cjoio- 1 * 

7.1. Raw materials and coal combustion plants 

The types of biomass waste used in this case study are FWR 
and AWR from Gomez et al. (2010a,b). The amount of biomass 
available is estimated by a hierarchy of potential approach 
that integrates physical, geographical and technical limita¬ 
tions, providing an upper bound for the potential estimation. 
Gomez et al. (2010a) use geo-referenced information to esti¬ 
mate the biomass production. The starting point of the present 
work are the UTM coordinates of selected collection areas and 
their biomass production, determined in Gomez et al. (2010a). 
The size of the areas are decided to provide enough biomass to 
allow a minimum installed capacity of hypothetical combus¬ 
tion or gasification plants at small scale. The reference areas 
are of 60 km x 60 km for AWR and 80 km x 80 km for FWR, 
regarding their territorial distribution and LHV. 

We assume that no transportation cost is charged for 
biomass assembly inside the collection areas. Only those areas 
that produce more than 50 kt compiled are considered for the 
co-combustion network. The amount of available biomass is 
reduced to take into account a certain competition with other 
matter applications. Fig. 3 depicts the location and the rela¬ 
tive size, in energy terms, of the selected biomass providers. 
Table 1 sums up the main characteristics of FWR and AWR 
used for modeling purposes. 

The selected technique that allows biomass usage in a com¬ 
bustion plant is co-firing. According to Basu et al. (2011), the 
investment is 192€/kW th . Operation and maintenance (O&M) 
costs are 4% the investment (Gomez et al., 2010b). The existing 
group of coal combustion power plants in Spain is considered 
by its used capacity during the year 2010, with an average 
working hours of 2800 h/yr (Red Electrica de Espana, 2010). 
See in Fig. 3 the power plants represented by their maximum 
percentage of inlet thermal power that can be replaced by 
biomass, i.e. 15%, without any efficiency penalty. The type of 
coal used into each plant is specified in Lopez-Vilarino et al. 
(2003) . See in Appendix A Tables 7 and 6 for further detail about 
the power plants, their coal origin and power produced. 

7.2. Seasonality and storage 

The total amount of energy contained in the biomass waste 
that can be used in the co-combustion power plants is set in 
44 PJ/yr; 15 PJ/yr coming from FWR and 29 PJ/yr from AWR. It is 


1 Chemical Engineering Plant Cost Index (CEPCI) from Chemical 

Engineering. 


assumed that the amount of coal that has to be replaced is a 
10% of the total inlet thermal power, 22 PJ/yr. The yearly energy 
profile is represented in Fig. 4: the period starts in March in 
order to allow biomass storage for the upcoming months. It is 
important to remark that between the biomass potential and 
the thermal demand, the SC efficiency, i.e. matter and energy 
losses in storage and pre-treatment echelons, “decreases” the 
initial available energy that have to match with the power 
plant requirement (see Fig. 5). 

Biomass storage is allowed after harvesting/collection and 
after pre-treatment sites, according to the general formula¬ 
tion. However, as pointed out in the previous section, this 
case study only considers storage or pre-treated matter. The 
present case study considers open air covered storage for 
raw material. MC decreases due to natural drying and DM is 




• • 


UTM coordinate x (zone 30 north) (km) 


Fig. 3 - Feedstock sites and the power plants location map. 
They are represented according to their relative GJ/yr 
production for sources, and according to their relative GJ/yr 
thermal input demand for power plants. Potential sites for 
pre-treatment and storage selection are selected among 
feedstock sites that produce higher amounts of biomass. 



7 8 9 10 11 12 1 2 

Month 


Fig. 4 - Biomass seasonal availability and power plants 
demand, according to a 10% of inlet thermal energy 
replacement. 
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reduced due to degradation: losses of 2% MC and 0.25% DM 
can be accounted (Rentizelas et al., 2009a; Maciejewska et al., 
2006). Raw biomass waste can be stored a maximum of two 
months. Storage of raw biomass allows for mixtures of the 
same type of biomass, gathering it from diverse sites and at 
different months. Mixtures of FWR and AWR are not allowed. 
It is assumed that storage of pre-treated biomass, in silos or 
tanks, does not change biomass properties. Therefore, this 
echelon, that preceeds power plants, is modeled as a flexible 
activity. 

7.3. Pre-treatment units and transportation 

Fig. 5 depicts the general network considered. After being col¬ 
lected, the biomass waste may be transported to different 
sites to be stored. Following, there are two mandatory pro¬ 
cesses before biomass pre-treatment to obtain the mandatory 
conditions of MC and shape: chipping and drying. Thereafter, 
torrefied biomass (TOR), torrefied pellets (TOP), pellets (PEL), 
bio-oil (OIL) or bioslurry (SLU) are produced before being stored 
if needed, and processed in the power plant. Trucks, adapted to 
carry solids or liquids are in charge of biomass transportation. 
The average working hours are 2400 h/yr. 

• Chipping. Based on a roll crusher, this unit increases the BD 
of raw biomass through conversion into chips to 240 k^m 3 . 
A loss of 2% of raw matter is homogeneously distributed 
among MC and DM (Hamelinck et al., 2003). The capacity 
ranges 10 and 80 ton/h. 

• Drying. A rotatory drum dryer is used to decrease the MC of 
the inlet mixture to 15%, 10% or 5%, depending on the pre¬ 
treatment unit requirement. Chipped biomass (from the 
previous echelon) changes its LHV and DM content accord¬ 
ingly (Hamelinck et al., 2005). The capacity ranges 40 and 
100 ton/h. 

• Torrefaction. Biomass with 15% MC is feeded to this unit. 
Torrefied biomass represents 70% of its initial weight and 
90% of its early energy content on a daf (dry and ash free) 
basis. The remaining humidity is 6% the inlet amount. The 


final BD is 230kg/m 3 (Uslu et al., 2008). The capacity ranges 
10 and 160 MW of the inlet energy content. 

• Pelletization. This activity processes biomass with 10% MC. 
The final bulk density is 575 kg/m 3 . The remaining humid¬ 
ity is a 17% of the inlet value. The energy loss is assumed 
to be 94% (Uslu et al., 2008). The capacity to be installed is 
between 10 and 160 MW of the inlet energy content. 

• Pelletization of torrefied biomass. TOP process increases the 
BD to 800kg/m 3 . MC decreases to 5% of the inlet biomass 
humidity (Uslu et al., 2008). The capacity is within 10 and 
160 MW of inlet energy content. 

• Fast pyrolysis. The unit processes biomass with a maximum 
5% MC. Bio-oil represents 73% of the initial weight and 66% 
of the inlet energy content on a daf basis. The remaining 
humidity is 64% the inlet value; the yield of water produced 
is between 5% and 15% on a dry basis. The BD results in 
1200kg/m 3 (Magalhaes et al., 2009; Uslu et al., 2008). The 
capacity is within 20 and 160 MW of inlet energy content. 

• Fast pyrolysis and char grinding. The bioslurry is 74% of the 
initial weight and 92% of the inlet energy content on a daf 
basis. The final MC is 53% of the inlet humidity. The BD is 
1200 kg/m 3 (Magalhaes et al., 2009; Uslu et al., 2008). This 
unit capacity ranges between 1 and 20 ton/h. 

• Transportation. Raw and pre-treated biomass is moved by 
trucks. Trucks for solid biomass can handle a maximum of 
130m 3 , while trucks for liquid biomass, 33m 3 . Moreover, 
trucks that move pyrolysis oil are more expensive than nor¬ 
mal carbon steel tanks (Hamelinck et al., 2003). Distances 
(in km) are assumed lineal from site to site with a tortuosity 
factor of 1.4 (Gomez et al., 2010b; Yu et al., 2009). 

See in Fig. 5, as example, the changes in properties of a 
stream of FWR. The network with the highest efficiency corre¬ 
sponds to pelletization of torrefied biomass, while the lowest 
value is for fast pyrolysis. 

The total capital requirement is spent at the beginning 
of the project. The scale factor is 0.7 (Uslu et al., 2008; 
Yu et al., 2009). Table 8 (see Appendix A) compiles invest¬ 
ment, O&M costs and utilities consumption for pre-treatment 
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Fig. 6 - Objective functions results for each SCN and 
extreme of the PF. 


units, storage and transportation echelons. Utilities price 
are 1393€/t (Ministerio de Industria, 2010) for diesel and 
0.04463 €/kWh (Comision Nacional de Energia, 2010) for elec¬ 
tricity. The electricity sold is bought at 0.03701 €/kWh (OMEL, 
2011). The depreciation period is 10 years: on the one hand, 
the combustion power plants have been constructed before 
2008, and on the other, the minimum average life for a pre¬ 
treatment unit is 10 years (Red Electrica de Espana, 2010). The 
candidate sites to allocate pre-treatment units and storage are 
selected among the locations with higher amounts of waste 
produced, above 95kt/yr, to avoid transportation costs (see 
Fig. 3). 

8. Results 

The following results show the capability of the developed 
model to propose optimal bio-based SC according to the 
economic and environmental criteria defined in Section 5. Eco¬ 
nomic optimization (EcO) and environmental optimization 
(EnvO) provide accordingly the results that correspond to the 
two ends of the Pareto frontier. The results cover: (i) objec¬ 
tive function values, (ii) breakdown of monetary spends, (iii) 
breakdown of environmental impact, (iv) technology installed, 
dimensions and flowrates among SC sites and echelons 
by month, and (v) elucidation about motivating policies or 
actions. Three scenarios are proposed to evaluate the trade¬ 
offs among the state-of-the-art (SOTA) pre-treatments, i.e. 
torrefaction, pelletization, pelletization of torrefied biomass, 
fast pyrolysis and fast pyrolysis combined with char grinding. 
The scenarios have been proposed to discern the best pre¬ 
treatment option, the best option for the case study (including 
the no use of any pre-treatment) and the best liquid-producer 
pre-treatment: 

• Scenario 1 (SCN1). This option analyzes all the SOTA 
pre-treatments proposed as option before biomass co¬ 
combustion, looking for the optimal technology selection. 

• Scenario 2 (SCN2). This scenario considers the optimization 
of the whole superstructure, including the no use of SOTA 
pre-treatments (however, chipping and drying are manda¬ 
tory). 

• Scenario 3 (SCN3). This alternative optimizes the SC consid¬ 
ering as SOTA pre-treatments only fast pyrolysis and fast 
pyrolysis with char grinding. 

8.1. EcO and EnvO 

Fig. 6 depicts the values of the two objective functions by 
scenario. The selected pre-treatments are, (i) torrefaction and 


pelletization of torrefied biomass for the economic criterion, 
only torrefaction for the environmental one, in SCN1, (ii) the 
no use of any SOTA pre-treatment for SCN2, and (iii) the use of 
fast pyrolysis combined with char grinding in SCN3 for both 
criteria. The NPV losses ranges between 0.7 and 2€ billions, 
while the environmental impact saving oscillates between 
270 and 425 thousands of pts/yr. At a first sight, SCN2 is 
the scenario that provides less NPV losses and higher envi¬ 
ronmental impact savings, followed by SCN1 and SCN3. It is 
remarkable that all the solutions derive into an environmen¬ 
tal impact reduction if compared to the use of only coal, even 
with the production of liquid biomass, which uses small trucks 
for product transportation, hence with higher environmental 
impact than the trucks used for solid transportation (see as 
reference Table 8 in Appendix A). As the biomass providers 
are highly dispersed along the territory in this case study, the 
trade-off between the investment in a pre-treatment technol¬ 
ogy vs. savings in transportation (due to the higher BD of the 
pre-treated biomass) favors the simplest choice. To verify and 
discuss this point, we analyze Figs. 7 and 8. 

In Fig. 7 the breakdown of initial investments and annual 
operation costs are depicted. From Fig. 7(a) it is recognized that 
for SCN1 and SCN3, almost half of the investment comes from 
the adaptation of combustion plants to fire biomass, while 
in SCN2 this is the main item. For SCN1 and SCN3 chipping 
and drying contribute to the investment less than the half of 
the costs of SOTA pre-treatments. Note that in environmental 
optimum solutions the investment in storage is higher than 
in economic ones. Fig. 7(b) shows that the difference between 
scenarios due to operational costs is less than the difference 
perceived for the investment. The major contribution comes 
from the acquisition of raw material. However, this item is 
slightly smaller for SNC2; the difference comes from the effi¬ 
ciency of the chains (see Table 9 in Appendix B), that takes 
into account the loss of matter and energy along the differ¬ 
ence processes from the collection step up to the entrance 
to the co-combustion plant. Differences among scenarios 
are coming from the installed technologies: a summation of 
utilities, O&M and transport. Concerning transportation differ¬ 
ences and storage investment differences observed previously 
between economic and environmental optimum configura¬ 
tions, it is deduced that the environmental optimum reduce 
impact by reducing transportation, and therefore, increasing 
storage. 

Fig. 8 depicts the environmental impacts resulting from the 
different bio-based SC activities. Fig. 8(a) reveals that the elec¬ 
tricity generation activity by biomass combustion is the item 
that contributes the most to the overall impact. This is fol¬ 
lowed by SOTA pre-treatment options. The lowest contribution 
comes from chipping, followed by transportation and biomass 
collection. According to Fig. 8(b), the main changes between 
economic and environmental optimum configurations are in 
transportation impact, followed by minimum differences in 
pre-treatment for SCN1, in transportation and drying for SCN2 
and just in transportation for SCN3. In SCN1 the differences 
in pre-treatment impact are due to the use of torrefaction and 
TOP pellets, or just torrefaction. In SCN2, the decrease in trans¬ 
portation impact is directly related to an increase of dryers 
installation. The contribution of drying and pre-treatments 
to the environmental impact in SCN3 (use of fast pyrolysis 
combined with char grinding) is larger than SCN1. 

Overall, SCN2 is more advantageous in investment and 
environmental impact, while the difference among scenarios 
is less important just regarding O&M costs. As suspected at 
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(a) Investment needed at the beginning of the project. 


EcO-SCNl EnvO-SCNl EcO-SCN2 EnvO-SCN2 

(b) Annual costs share. 


Fig. 7 - Breakdown of costs for economic and environmental optimal networks for each SCN. 



(a) Environmental impact distribution among all the 
SC activities. 

Fig. 8 - Distribution of the environmental impact by echelon 
environmental optimal networks. 



EcO-SCNl EnvO-SCNl EcO-SCN2 EnvO-SCN2 EcO-SCN3 EnvO-SCN3 


(b) Share of the environmental impact distribution 
without considering electricity generation. 

1 bio-based SC for each SCN and economic and 


the beginning of the section, for this specific case study, the 
extra-investment in SOTA pre-treatments is not justified by 
the environmental impact of transportation. 

See in Table 2 the comparison, in terms of environmen¬ 
tal impact, between the contribution of the biomass SC and 
the coal replaced SC for each optimum configuration and sce¬ 
nario. It is seen that the coal impact is different even if inlet 
energy replaced is the same in all the cases. This is due to the 
fact that the inlet coal replaced per power plant is different 
into each configuration (see next subsection). The most con¬ 
tributing impact categories are climate change for biomass, 
and human health followed by climate change for coal. Note 
that the replaced coal SC’s that have higher impacts corre¬ 
spond to environmental optimizations. Accordingly, biomass 


SC’s that optimize the economic criteria have higher impact 
than environmental optimum SC’s. 

8.2. Optimum network’s description 

Figs. 9-11 illustrate the optimum networks obtained for each 
criterion and scenario. Biomass providers are characterized 
by their relative GJ/yr used from each production site to sup¬ 
ply the demand. Co-combustion plants are represented by 
their replaced thermal input in GJ/yr. The lines reproduce 
transportation of biomass between sites and echelons, implic¬ 
itly, in terms of ktons/month. It distinguishes transportation 
of raw matter from suppliers toward intermediate (pre¬ 
treatment) sites, transportation of pre-treated matter between 


Table 2 - End-point impact categories for each optimum c< 
SC corresponding to the replaced coal. 

mfiguration, taking 

into account the 

only biomass SC and the 

End-point impact category 

EcO-SCNl 

EnvO-SCNl 

EcO-SCN2 

EnvO-SCN2 

EcO-SCN3 

EnvO-SCN3 

Biomass 

Human health 

56,202.98 

52,397.64 

27,350.03 

27,095.71 

58,077.93 

56,683.34 

Ecosystem quality 

9027.57 

8294.78 

6815.16 

6591.74 

9927.48 

9613.97 

Climate change 

421,785.87 

419,046.10 

365,294.84 

357,574.40 

442,714.77 

441,636.23 

Resources 

32,233.37 

28,636.16 

17,254.91 

19,029.81 

37,123.09 

35,931.63 

Impact 2002+/yr 

519,249.79 

508,374.69 

416,714.94 

410,291.66 

547,843.27 

543,865.17 

Coal replaced 

Human health 

379,305.62 

385,769.57 

383,845.49 

385,774.88 

383,567.11 

385,771.91 

Ecosystem quality 

20,438.46 

22,539.02 

21,914.52 

22,541.43 

21,824.01 

22,540.08 

Climate change 

232,909.63 

240,757.59 

238,417.47 

240,760.41 

238,079.73 

240,758.83 

Resources 

176,593.30 

185,578.79 

182,928.67 

185,607.99 

182,540.19 

185,591.66 

Impact 2002+/yr 

809,247.01 

834,644.97 

827,106.16 

834,684.71 

826,011.04 

834,662.48 
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I Table 3 - Total capacity installed in the optimum configurations. I 


EcO-SCNl 

EnvO-SCNl 

EcO-SCN2 

EnvO-SCN2 

EcO-SCN3 

EnvO-SCN3 

Chipper (ton/h) 

974.07 

1029.89 

903.07 

1023.34 

1029.78 

1115.00 

Dryer (ton/h) 

952.26 

1906.14 

880.47 

990.00 

1137.62 

1180.00 

Torrefaction (MW t hi„) 

1070.00 

3062.12 





Pelletization of torrefied 
biomass (MW t hi n ) 

Fast pyrolysis and char grinding 

1875.18 




681.22 

715.83 

(ton/h) 

Storage matter 

Torrefied pellets 

Torrefied biomass 

Dried chips (15%) 

Dried chips (5%) 

Slurry 

Slurry 

Storage (m 3 ) 

19,005.41 

153,142.37 

27,821.21 

164,971.09 

6208.06 

10,110.84 


intermediate sites (toward SOTA pre-treatment or storage) 
and transportation of pre-treated matter to the co-combustion 
plants. For a connection to appear in the map, it is enough to 
be presented just a month. It means that the transportation 
depicted refers to movements by year, occurring or not at the 
same month (see Figs. 13 and 14 in Appendix B.l for clarifica¬ 
tion). Even if the connection among sites are represented by 
straight lines, the case study is considering a tortuosity factor. 
Table 3 summarizes the total capacity installed per technology 
for each optimum network. 


Regarding Fig. 9(a) and (b) the most characteristic feature 
is the difference among the type of transported matter: in 
the second one no transportation between intermediate sites 
is appreciated. Moreover, the covered co-combustion plants 
are clearly different. As observed in the previous section, the 
environmental optimum scenario reduces the environmen¬ 
tal impact by lowering transportation or connection between 
nodes. Considering the capacities installed as in Table 3, 
chipping, drying and SOTA pre-treatment capacities are very 
similar while the capacity installed for storage is around 



(b) Optimum Impact 2002+ network configuration for SCN1. (b) Optimum Impact 2002+ network configuration for SCN2. 


Fig. 9 - Optimum network configurations for the two 
selected criteria and mandatory SOTA pre-treatments 
(SCN1). 


Fig. 10 - Optimum network configurations for the two 
selected criteria and only chipping and drying as 
pre-treatments (SCN2). 
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(a) Optimum NPV network configuration for SCN3. 



UTM coordinate x (zone 30 north) (km) 

(b) Optimum Impact 2002+ network configuration for SCN3. 


Fig. 11 - Optimum network configurations for the two 
selected criteria and production of liquid fuel (SCN3). 

8 times higher in the environmental optimum network. See 
also in Table 3 the type of biomass that is stored. The 
results from the model explicites that the number of sites 
chosen to store biomass for the environmental optimum net¬ 
work is higher and with larger capacities per site than for 
the economic optimum network. In contrast, in the eco¬ 
nomic optimum network, there are less pre-treatment units 
installed, at their maximum capacity, while for the other opti¬ 
mum, in almost every candidate site there is a dryer installed 
at its minimum available capacity. According to Figs. 6 and 9, 
it is cheaper to transport less matter and install more concen¬ 
trated capacities, since economies of scale are applied. Storage 
has a certain investment cost, but no O&M costs neither envi¬ 
ronmental impact: the environmental optimum network does 
not limit the installation of storage sites. 

For Fig. 10(a) and (b) about SCN2, analogous features can 
be remarked as in SCN1 by regarding transportation and 
technology installation dispersion. However, the covered co¬ 
combustion plants are more homogeneous in SCN2 than in 
SCN1 optimum networks. Note that in this scenario the stor¬ 
age capacity installed is higher than the capacity installed in 
SCN1: chipped biomass is less dense than torrefied and pel¬ 
letized biomass, so storage is larger even if less biomass is 


used due to the SC efficiency (see Tables 3 and 9 (see Appendix 
B)). Note that slightly more capacity of chipping and drying 
is installed when the SC efficiency is lower, as remarked in 
Table 9 (see Appendix B). 

Fig. 11(a) and (b) depict the optimum networks for fast 
pyrolysis combined with char grinding. Both networks are 
very similar, with the exception of specific connections and 
co-combustion plants coverage. There is a predominant trans¬ 
portation of raw matter from supplier to intermediate, even 
if biochar is more dense: transportation by trucks for liq¬ 
uid biomass is done by smaller units. The storage capacity 
installed in this scenario is significantly less than in the other 
ones. From the size of the fast pyrolysis plants installed and 
also contrasting with the efficiencies of chipping and drying 
(as for example calculated in Fig. 5), it can be deduced that 
the capacities installed for chipping and drying (for the three 
scenarios) are larger than needed. It puts into manifest a “mis¬ 
match” between the needed capacities and the assumed exist¬ 
ing capacities in the market; therefore, there is installed capac¬ 
ity that is not used. On the other hand, fast pyrolysis units are 
installed in around 70% of the candidate nodes at its maxi¬ 
mum capacity: it means that the granularity of the solution 
is probably higher than needed due to another “mismatch” 
between the needed capacity and the available capacities. 

According to the biomass distribution along the territory 
in Fig. 3 and to the detailed amounts of biomass used in 
Table 9 (see Appendix B), the optimal solutions use first the 
biomass with higher LHV, that is also the non-seasonal type, 
and then, the closer available spots. The model uses all FWR 
and adapts the consumption of AWR to cover the demand of 
the co-combustion plants regarding the SC efficiency. For that 
reason, SCN3 is consuming more biomass. 

Fig. 12 characterizes the portion of coal that is replaced by 
biomass into each plant, distinguishing among scenario and 
optimum network. Plants m7 and ml3 do not produce energy 
according to the selected input data. Plants mil, ml6 and ml7 
use imported coal. It can be recognized that the model selects 
the inlet energy to cover by plant independently of the size of 
the plant (see Fig. 3 and/or Table 7, see Appendix A). However, 
it considers the type of coal consumed, i.e. local or imported, 
and plant location. The selected plants in all the optimum 
networks to replace their coal at the maximum value (15%) 
are ml2, ml4, ml5, ml8 and ml9. ml, m2 and ml6 have also 
coal replaced by biomass in all the optimum networks, but 
at a variable portion depending on the scenario, m3 and m8 
are the less used plants. Plants that use imported coal have it 
replaced by biomass at the maximum fraction in EnvO chains, 
ml, which is the plant that produces the largest amount of 
energy, is always selected to be replaced at 15% by the EcO 
chains. SCN2 optimum configurations use less plants than the 
other two scenarios. 

Figs. 13 and 14 (see Appendix B.l) represent the informa¬ 
tion, i.e. mass flowrates by month of raw biomass, pre-treated 
biomass that goes to storage or pre-treated biomass that goes 
to a co-combustion plant, that can be retrieved for each active 
node in the optimum solutions, using two specific sites as 
example. 

8.3. Discussion 

The previous results shows that: 

• The difference between EcO and EnvO networks comes from 
the installation of more or less dispersed pre-treatment 
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■ EcO-SCNl ■ EnvO-SCNl ■ EcO-SCN2 ■ EnvO-SCN2 ■ EcO-SCN3 ■ EnvO-SCN3 



Fig. 12 - Portion of electricity produced by coal (maximum of 0.15) that is replaced by biomass. 


units (less advantage of economies of scale) and the use 
of transportation. 

• The scenario with liquid biomass points out a use of sig¬ 
nificantly less storage. Moreover, the costs associated to 
transportation in SCN3 are higher than in SCN1, but lower 
than in SCN2. 

• Results are highly influenced by the available unit’s capacity 
and the costs associated to each echelon. 

• This specific case study has a well distributed raw matter 
along the territory. The model chooses to use the residue 
with higher LHV and without seasonality. Large transporta¬ 
tion distances are not needed. 

The previous results have also elucidated important trade¬ 
offs to be analyzed when aiming to implement a biomass SC in 
any case study, i.e. the selection of the SOTA pre-treatment(s): 

• Transportation costs and environmental impact vs. 
installed capacities. 

• Granularity vs. capacities existing in the market. 

• Storage size vs. transportation costs. 

• Distances to cover vs. transportation costs and environmen¬ 
tal impact. 

The following paragraphs analyze the inlet data and the 
results obtained for SCN1, SCN2 and SCN3 in order to favor the 
use of SOTA pre-treatments, to (i) motivate the conversion into 
liquid biomass and to (ii) persuade the use of solid enhanced 
biomass, such as torrefied biomass and pellets. The last anal¬ 
ysis aims at calculating the electricity price that would allow 
to invest in such a project. 

8.3.1. Enhanced biomass 

Due to its higher BD, liquid biomass is adequate for long 
distance transportation and for storage in smaller places. Nev¬ 
ertheless, fast pyrolysis is still in a pre-commercial state and 
therefore, it is much more expensive than other alternatives. 
In addition to that, according to the data for this specific case 
study, trucks to move liquid biomass are smaller and more 
expensive than trucks used to transport solid biomass. If com¬ 
paring just transportation costs associated to EcO and EnvO 
of SCN1 and SCN2, it is seen that this item should decrease 


around 30-40% to equiparate the different technologies, or 
from another point of view, the distances to be covered should 
increase in a 1.4-1.5 to justify the conversion into biochar. If 
comparing (i) consumed biomass (i.e. efficiency of the pro¬ 
cess), (ii) transportation and (iii) investment costs for the same 
networks, it results in a needed global cost decrease of 5-12% 
to promote the biochar; specifically with a 15-28% of invest¬ 
ment decrease. Another possibility to motivate liquid biomass 
use is to restrict the amount of installed storage by a legislative 
tool. 

If comparing SCN3 and SCN2 for liquid biomass promotion, 
the results are more extreme. The investment in technology 
should decrease in around a 60%. Moreover, a change just in 
transportation costs or distances, will not favor the use of liq¬ 
uid biomass if compared to the use of no SOTA pre-treatments. 
On the other hand, looking at the environmental impact, 
the difference between transportation impacts for SCN3 and 
SCN2 puts into manifest a needed increase of 5.5-7.6 to jus¬ 
tify biomass conversion into biochar. This can be motivated 
by changing the means of transportation, i.e. using pipeline, 
train or ship, and therefore, implying longer distances. 

By comparing SCN1 and SCN2, the investment in the 
pre-treatment technology should decrease around 30-45% to 
motivate the use of torrefied biomass or TOP pellets. In an 
analogous way as for the previous comparison, the transporta¬ 
tion impact should decrease 5.5-7 times to favor the use of 
enhanced solid biomass. 

8.3.2. Price 0/kWh 

See in Table 4 the needed increment in the sold electricity 
price to afford the investment and operation along 10 years 
of each one of the optimum networks proposed for each sce¬ 
nario; it oscillates between 9 and 23%. This monetary charge 
can be relieved by specific subsidies like grants, tax incen¬ 
tives, bonuses, feed-in tariffs or emissions trading, according 
to specific policies to decrease the environmental impact in 
the electricity generation sector. 

9. Conclusions 

This paper presents an approach for efficient design and 
planning of bio-based SC. The model relies on multi-period 
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EcO-SCNl EnvO-SCNl EcO-SCN2 EnvO-SCN2 EcO-SCN3 EnvO-SCN3 

Revenues increase 15.30% 20.66% 9.20% 11.08% 21.71% 23.24% 


MILP for multi-objective optimization of economic and envi¬ 
ronmental interventions. The model considers long-term 
strategic decisions, such as the selection of biomass sources, 
establishment of pre-treatment units and their location and 
disposition of distribution centers. The specific case study 
contemplates a bio-based SC located in Spain and includes 
an estimation of the potential of woody residues (FWR and 
AWR) to supply the coal power plants that exist in the country. 
The main purpose in this context was to provide the opti¬ 
mum SC networks in the light of different scenarios that affect 
the pre-treatment choice to replace a 10% of the coal inlet 
thermal power. Economic and environmental criteria have 
enabled the discrimination of the optimum networks, allow¬ 
ing specific guidelines for the orientation of economic and 
environmental policies. It has been demonstrated that for the 
specific case studies and environmental metric considered 
in this work, co-combustion is always positive to decrease 
environmental impact. These results are based on the under¬ 
lying assumptions of the LCA boundaries. The environmental 
impact reduction depends on the specific raw material com¬ 
position and on the specific environmental effects considered, 
which for the case of the cofiring case are assumed to be the 
same as those of coal. 

Even though the SOTA pre-treatment processes proposed 
in this paper have been developed, used, and are contin¬ 
uously improved, its use nowadays influences (penalizes) 
the SC efficiency and cost. However, their implementation 
enhances biomass properties as well as homogenizes the raw 
matter. An important advantage would be the use of pre¬ 
treated biomass to cover large distances (international SC) or 
to move large amounts of biomass. According to the results, 
in a national or local scale with sources of biomass waste 
well distributed along the territory, the installation of pre¬ 
treatments would obey a policy, strategic decision to dispose 
organic matter or decrease the environmental impact. SOTA 
pre-treatments will have to evolve toward more efficient tech¬ 
nologies, more integrated processes with lower investment 
and O&M costs. It has been seen that the model encourages the 
use of the biomass that has a higher LHV beyond biomass that 
is closer to the destination point. Overall, the selection of the 


simplest SC is recommended when large distances do not 
need to be covered or when there are no storage space restric¬ 
tions. In that case, liquid biomass would be encouraged. 
Storage of pre-treated biomass, liquid or enhanced solid, is 
adviced for longer periods of storage that may affect the prop¬ 
erties of the biomass stream, i.e. matter degradation, MC 
acquisition or even microorganisms formation. The use of a 
SOTA pre-treatment can be also motivated by biomass final 
process specific requirements. 

This work has solved a specific case study and the results, 
i.e. the proposed optimum networks, the specific environ¬ 
mental impact and economic values, the changes that have 
to occur to motivate different technologies and the needed 
revenue increase to support such a project, cannot be extrap¬ 
olated. However, the same approach and analysis can be 
performed for other input data, other objective functions, and 
other case studies. SC are tailor-made and are influenced by 
the criteria of the decision-maker. 

Future work will have to focus on more exhaustive biomass 
properties calculation, including biomass different types 
combinations, long-term storage and seasonality. A more 
exhaustive electricity demand can be estimated, taking into 
account the fluctuation among years, not only seasonality, as 
well as uncertainty. Calculation time can be a limitation in a 
bio-based SC problem resolution. Decomposition methods are 
envisaged to tackle with it. The characterization of the Pareto 
Fronts for each scenario is also proposed. 
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Table 6 - Types of coal used in the Spanish power plants. Data for 2010 (Club Espanol de la Energfa, 2010; Lopez-Vilarino 
et al., 2003; Red Electrica de Espana, 2010). 


Type 

Origin 

LHV„ r (MJ/kg) 

Cost (€/t) 

Coal 1 

Sub-bituminous 1 

Local 

12.57 

85 

Coal 2 

Sub-bituminous 2 

Local 

17.81 

85 

Coal 3 

Bituminous 1 

Local 

22.63 

85 

Coal 4 

Bituminous 2 

Imported 

27.03 

80 


Table 7 - Spanish thermal power plants characteristics ordered by region, 
for 2010 (Lopez-Vilarino et al., 2003; Red Electrica de Espana, 2010). 

Installed power (MW) and used capacity (GWh) 

Model name 

Name 

Power (MW) 

Energy (GWh) 

Type of coal 

ml 

Puentes Garcia Rodriguez 

1468 

4955 

Coal 1 

m2 

Meirama 

563 

856 

Coal 1 

m3 

Abono 

916 

3663 

Coal 3 

m4 

Lada 

513 

698 

Coal 3 

m5 

Soto de la Ribera 

604 

927 

Coal 3 

m6 

Narcea 

595 

1 

Coal 3 

m7 

Anllares 

365 

0 

Coal 3 

m8 

Compostilla 

1171 

209 

Coal 3 

m9 

La Robla 

655 

29 

Coal 3 

mlO 

Guardo 

516 

63 

Coal 3 

mil 

Pasajes de San Juan 

217 

487 

Coal 4 

ml2 

Ceres 

162 

516 

Coal 2 

ml3 

Escatron 

80 

0 

Coal 2 

ml4 

Teruel 

1102 

1793 

Coal 2 

ml5 

Escucha 

159 

156 

Coal 2 

ml6 

Litoral de Almeria 

1159 

4409 

Coal 4 

ml7 

Los Barrios 

589 

2489 

Coal 4 

ml8 

Puertollano 

221 

255 

Coal 3 

ml9 

Puentenuevo 

324 

590 

Coal 3 

Total 


11,379 

22,096 



Table 8 - Economic parameters of pre-treatment and storage units and transportation (Hamelinck et al., 2003; Magalhaes 
et al., 2009; Uslu et al., 2008). 


Base scale 

Base investment 

O&M (% of investment) 

Utility 

consumption 

Lifetime 

(yr) 

Chipping 

80t/h 

0.7 M€ 

20 

Bond law 

15 





0-15 • tj„ plJ t (kW) 


Drying 

100 t/h 

6.9 M€ 

3 

0.06 • tH 2 oeu (t diesel) 

15 

Torrefaction 

40MW thjn 

6.2 M€ 

5 

92 kWh/ti n p Ut 

10 

Torrefaction + pelletization 

40MW thjn 

7.5 M€ 

5 

102 kWh/ti n p Ut 

10 

Pelletization 

40MW th ^ 

5.9 M€ 

5 

129 kWh/ti n p Ut 

10 

Fast pyrolysis 

40MW thjn 

10.5 M€ 

4 

75 kWh/ti n p Ut 

25 

Fast pyrolysis + char grinding 

5t/h 

4.9 M€ 

5 

95 kWh/tj n j, ut 

25 

Open air covered storage 



0.53 €/m 3 month 



Silo storage 

5000 m 3 

0.45 M€ 

3 


25 

Tank storage 

2272 m 3 

1.11 M€ 

3 


25 

Transportation 

130 m 3 


0.69 €/m 3 



Solid 



1.16€/km 



Transportation 

33 m 3 


0.69 €/m 3 



Liquid 



1.70 €/km 
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Appendix B. Output data from the resolution of 
the case study 


Table 9 - SC efficiency, defined as the calorific value between biomass that is transported to the co-combustion plant, and 
raw matter. Woody residues consumed by optimum network and scenario. The maximum production of each residue is 
1162.20 ton for FWR and 2717.80 ton for AWR. 


EcO-SCNl 

EnvO-SCNl 

EcO-SCN2 

EnvO-SCN2 

ECO-SCN3 

1 EnvO-SCN3 

Eff (%) 

85.53% 

85.38% 

93.68% 

95.34% 

81.76% 

81.76% 

FWR consumed (kton) 

1162.20 

1162.20 

1162.20 

1162.20 

1162.20 

1162.20 

AWR consumed (kton) 

1046.94 

1051.20 

838.75 

800.85 

1157.13 

1157.13 


B.l. Information available for each site 

The following figures represent the information that can be 
retrieved for each biomass provider and/or pre-treatment site 
that is active in the optimal solutions, i.e. the mass balance 
by month. This information includes raw biomass compiled, 


production of pre-treated biomass, inventory, compilation 
of pre-treated biomass, pre-treated biomass that is sent to 
another nodes, and pre-treated biomass that is sent to co¬ 
combustion plants. See as example the planning for nodes/cl 
and/c2 for SCN1 and EcO (Figs. 13 and 14). 



(a) Raw biomass (FWR) that /cl receives (b) Production of pre-treated biomass by /cl. 
from other sites. 



(c) Pre-treated biomass (chipped FWR in (d) Pre-treated biomass (dried FWR 
green and chipped AWR in brown) that 15%) that /cl sends to other sites. 

/cl receives from other sites. 



(e) Pre-treated biomass (torrefied 
biomass) that /cl sends to co-combustion 
plants. 

Fig. 13 - Mass balance of site/d, i.e. planning of biomass flows. 
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(a) Raw biomass (FWR) that /c2 receives (b) Production of pre-treated biomass by /c2. 
from other sites. 



Month Month 


(c) Pre-treated biomass (chipped FWR) (d) Storage of pre-treated biomass, 
that fc 2 receives from other sites. 



(e) Pre-treated biomass (pellets of tor¬ 
refied biomass) that /c2 sends to co¬ 
combustion plants. 


Fig. 14 - Mass balance of site/c2, i.e. planning of biomass flows. 
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